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We compute the four-point correlation function of a light 0{N) scalar field in de Sitter space in 
the large- A*' limit. For superhorizon momentum modes, infrared effects strongly enhance the size of 
loop contributions. We find that in the deep infrared limit, the latter are of the same order as the 
tree-level one. The resulting nonGaussian correlator presents a singular momentum dependence. 
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Quantum field theory in curved spaces is a topic of 
great interest with a long history The case of de 
Sitter (dS) space has attracted a lot of attention both 
because of its large degree of symmetry and because of 
its phenomenological relevance for the early inflationary 
era and for the current accelerated expansion of the uni- 
verse. Specific phenomena such as gravitational redshift 
or particle creation imposes one to rethink much of what 
is known in Minkowski space, starting from the basic no- 
tions of particle and vacuum state, even for free fields j^l- 
At present, free gauge fields, such as the photon or the 
graviton, are still the subjects of debates [3]. 

Interacting fields can be studied by means of pertur- 
bation theory 0-[8|. They pose practical and conceptual 
issues. An example is the transPlanckian problem 
i.e., the question of the effective decoupling between in- 
frared (IR) and ultraviolet physics, which underlies the 
very concept of QFT on dS. They also reveal novel spe- 
cific features as compared to the flat space case. For 
instance, scalar fields of sufficiently large mass — in units 
of the expansion rate — are fundamentally unstable and 
can decay to themselves Light fields, which have 

no Minkowski analog, are also of great interest because 
of their phenomenological relevance, e.g., for infiation- 
ary cosmology. They exhibit strong semi-classical fiuc- 
tuations for superhorizon modes and turn out to be es- 
sentially nonperturbative, even at weak coupling, due to 
large IR effects [1, [ll|. In recent years, various methods 
inspired from flat space techniques have been developed 
to deal with IR issues in dS. Results are still rather scarce 
but the nonperturbative aspects of light scalar fields are 
being unravelled [l2| - [20j 

A typical example is the phenomenon of dynamical 
mass generation: a field with vanishing tree-level mass 
develops an effective mass due to its self-interactions [ll| . 
This lifts the fiat tree-level potential and regulates pos- 
sible IR divergences. Incidentally, this results in nonan- 
alytic coupling dependences of physical observables. A 
similar phenomenon has been demonstrated for an 0{N) 
scalar field in the large-iV limit in the case where the 
tree-level potential shows spontaneous symmetry break- 
ing [l5[. Strong IR fluctuations restore the symmetry, 
as anticipated in [2l| . and lead to nonperturbatively en- 



hanced loop contributions fl9'|. 

Immediate phenomenological implications of nontrivial 
field interactions in the infiationary universe are possible 
quantum corrections to standard inflationary observables 
,6, 22^, or the possibility of nonGaussian features of pri- 
mordial density fluctuations p3[ . It has been pointed 
out in [l^ that IR effects may lead to parametrically en- 
hanced nonGaussianities at tree-level. In this letter, we 
show that loop corrections to nonGaussian correlators 
are also amplified by IR effects and may in fact be as 
large as the tree-level contribution. We consider the case 
of an 0{N) theory with quartic interactions and employ 
a nonperturbative large- limit. Using the expressions 
for the field propagator and four-point vertex function 
recently obtained in Refs. [l5l.[l9|. we compute the equal 
time four-point correlation function for the phenomeno- 
logically relevant case of superhorizon modes, which we 
obtain in closed analytical form. We point out the inter- 
esting momentum structure of the loop contributions. 

Consider the 0(iV)-symmetric scalar field theory with 
classical action (a sum over a = 1, . . . , iV is implied) 



(1) 



with the invariant measure = J d'^'^^x \/ —g, on the ex- 
panding Poincare patch of a d + 1-dimensional dS space. 
In terms of comoving spatial coordinates X and confor- 
mal time — oo < ?7 < 0, the line-element reads (we choose 
the Hubble scale H = V) 



ds^ = r]-^ {-drf + dX ■ dX) 



(2) 



In Eq. ([T]), the mass term m^g — m? + £JZ includes a 
possible coupling to the Ricci scalar TZ — d{d+ 1) and □ 
is the appropriate Laplace operator. 

In the following we consider the n-point correlation 
and vertex functions of the conformaly rescaled fields 
4>a{x) = (— ?/)~2~(^^(j;) in the (interacting) Bunch Davies 
vacuum state. The latter are conveniently expressed 
in terms of time ordered products of field operators 
along a closed contour in (conformal) time; see, e.g., 
p^ . For instance the two-point function Gab{x,x') = 



2 




FIG. 1: The self-energy in the limit N ^ oo; see Eq. (fTO)l . 
The internal line corresponds to the propagator G itself, hence 
the nonperturbative character of this limit. 



{Tc4>a{x)(l)b{x')) , where Tq denotes time-ordering along 
the contour C, encodes both the statistical and spec- 
tral correlators Fab{x,x') = ^{{(j)a{x),(l)i,{x')}) and 

Pabix,x') = i{[(l)a{x),(l)bix')]): 

Gab{x,x') Fah{x,x')-^signc{x° -x'°)pab{x,x') , (3) 

where the sign function is to be understood on the con- 
tour C. It was shown in that, in the large- iV limit, 
the system only admits 0(A'')-symmetric solutions. We 
thus have ((/>a) = and Gab — 5abG. 

In the symmetric phase, the four-point correlation and 
vertex functions G'^'*^ and F*^^^ are related by 



GaBCD ^ GAA'GBB'GcC'GDD'i^A^B'c'D' 



(4) 



where capital letter indices collectively denote space- 
time variables and 0{N) indices and an appropri- 
ate integral/summation over repeated indices is under- 
stood. Here, we are interested in computing the equal- 
time four-point correlator in comoving momentum space 
G^^^ (rj, Ki, ... , K4) for superhorizon physical momenta, 
—KiTj < 1. Both the propagator G and the vertex F'^''^ 
have been computed recently in the IR regime in the limit 
TV ^ 00 [H [la. Let us briefly review the results relevant 
for our present purposes. 

In comoving momentum space, the two-point function 
has the free-field-like expression, for sign^{ri — rj') — 1, 



(5) 



where Hi,{z) is the Hankel function of the first kind and 
1^ — ^ (P- — A'P. Here, M a self-consistent, dynam- 
ically generated mass, to be discussed shortly. In the 
cases of interest below, M <C 1 and it is convenient to 
introduce the small parameter e — d/2 — v w KP /d. 
For superhorizon modes, the statistical and spectral two- 
point function read 



(6) 



Pm{v,v',K) = -VWn"(ln^), (7) 
where F^, ~ [2''F(j/)]^/47r and we introduced the function 



V'aix) 



sinh(aa;) 



-b\x\ 



(8) 



The self-consistent mass M satisfies the gap equation 

(9) 



= mis + 



where the constant a is given by the tadpole diagram of 
Fig. [TJ Retaining only the dominant IR contribution in 
the loop (see [l^ for a complete treatment), one gets 



(10) 



where we introduced Aoft = \F^VI^/\2{2ttY and Vl^ 
2TT'^/'^/T{d/2). Equation © is solved as 



'(™ds)' 



dX 



cff ■ 



(11) 



The nonanalytic coupling dependence reflects the non- 
perturbative IR character of mass generation. 

The four-point vertex function can be written as [l^ 

riLfe,K,) = [77i---7y4]^ 

>^\5ab5cd5c{vi-m)6c{m~Vi)iD{rii,ri:i,Ki2)+Y>crm}^, 

(12) 

where (5c (?7 — ??') is a Dirac delta function on the contour, 
Kij = \K.i + Kjl and 'perm.' denotes the two permuta- 
tions needed to make F^"*) symmetric. The function D is 
the two-point correlator of the composite field x 0^ : 

iD{ij, 7j', K) = -^ [6c{ii - 11') + ^liv, V',K)] , (13) 

The first term on the right hand side corresponds, when 
inserted in Eq. , to the tree-level vertex and the func- 
tion / resums an infinite series of bubble loop diagrams. 
The one-loop contribution gives a nonlocal power in the 
IR. Each additional loop leads to an additional power of a 
large IR logarithm which actually resum to a modified IR 
power law J^] . The nonlocal function / admits a decom- 
position in terms of statistical and spectral components 
as in ([3]), which read, in the IR, 



F, 



In^ 



(14) 
(15) 



with 7r„ 



2(7, 



7r„ and k 



e. Here, 7r„ 



the effective parameter which actually controls the loop 
expansion. In the following we shall assume TTp < 1 and 
we introduce the parameter s = d/2 — D. Expanding 
the above expression in powers of the effective coupling 
TTp generates the whole series of perturbative logarithms. 
One sees however that the latter breaks down for large 
time separations TTp | In 77/77' | > 1 and/or deep IR mo- 
menta TTp I \nK^'q'q'\ > 1. 
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function / in (IT3l). We write 



+ perm. 



K2 K3 

6 c 



FIG. 2: The equal-time four-point correlator (|16|) . The black 
lines represent the propagator G and the dashed line rep- 
resents the nonlocal vertex function (|12|) . which includes the 
tree-level vertex and an infinite series of bubble loop diagrams. 



We now have all the ingredients for our computation 
of the contribution from superhorizon, IR modes to the 
four-point equal-time correlator ([5]). Writing = /c '^Ci 
the latter can be expressed as the following integral on 
the time contour C 



(16) 



-|-perm. 

where we introduced the functions 

A,,{t^,0 = G{r,,^,K,)G{r^,^,Kj), (17) 

The tree- level contribution is cx Sc{£, — £,') and must be 
performed separately. It involves the time integral 

Jc 

d^-O'-' {AUv,OA',,{v,0 + (12 ^ 34)} ,(19) 



where we introduced the statistical and spectral compo- 
nents of the function Aij , as in ^jL and we used standard 
manipulations on the contour [24| . Note the symmetry 
relations (77, C) = Af^ (C, rj) and A^^ (77, ^ = -^f, v)- 
Due to the strong IR enhancement of the statistical func- 
tion ^ as compared to the spectral one ([7|), we have 



(20) 



A^MO = F{^,tK,)p{f^,C,K,) + {t^j), (21) 

where we neglected a term cx pp in the first line. This is 
typical of the classical statistical field regime 0, HI] and 
reveals, in the present context, the classical stochastic 
nature of dS IR fluctuations [26| . To estimate the contri- 
bution from superhorizon modes we replace the integral 
ri _^ where nn is such that the relevant momenta 

J— 00 -'^70 

are superhorizon: —Kirjo < 1. One can then use the 
expressions ([6]) and ([7]) to compute (IT9t : see below. 
The loop contribution in (jl6p involves the nonlocal 



—00 

7] 



di'AU^,OiF{i,i')K,{e,v) 



/ d£,'A!^{ri,Olp{^,aKi{£.'.v) 



di / rf?'Af;(r,,e)/p(e,e')^?,(e',^), m) 

-00 



replace again 

^ p and use the IR behaviors (|14p 
and (fT5|) . The calculation is straightforward. 

Extracting a overall factor and introducing the variable 
X = ln(7;/77o), our final result reads 

^(4) , A F3 (_^)2-4.(_^^)2e 



3N 2v [Ki--- K^Y'^ 
-l-perm. , 

with the two momentum structures 



(23) 



g{x,Ki) = g^{x) {Kf ' + ■■■+ Kf) 

The function gi{x) receives contributions from the tree- 
level vertex and from the last two lines of Eq. \22\ , while 
g2{x) comes from the second line of Eq. (|22|) . We find 



gi[x) = £2e(a;) H 



211 £ — £ 

52(x) = g(-r;o)2^£?+,(a:), 



(25) 
(26) 



where we defined the function Ca{x) = (e"^ — l)/a. Note 
the mild dependence on r^o, which indicates that our as- 
sumption of IR dominance is consistent. The first term 
on the right hand side of Eq. ([25|) is the tree-level contri- 
bution. Loop terms are cx TTp. 

Nontrivial IR effects arise in the cases of vanishing, or 
negative tree-level square mass TTt^g < 0. For TT^^g = 0, 
the dynamically generated mass pT|) is cx a/A and 



£ = i/Acff/d , TTp = 2de , 

such that £ = 3£. We thus get 

gi{x) = C2e{x) {I + eC2e{x)) , 
g2{x) - 2e{-i^of'CUx). 



(27) 



(28) 
(29) 



We see that loop corrections get enhanced by IR effects 
and are controlled by £ '-^ \/A. In the regime e\x\ < 1, 



gi{x)^x + £x\ g2ix) 2e{~r]of''x'^ 



(30) 



and we recover the usual perturbative result [13]: G*-^-* ^ 
Xx{l+0{V^x'^)). However, we see that for A]a;] ~ 1, loop 



4 



contributions become of the same order as the tree-level 
one. In the deep IR regime e\x\ > 1, the linear growth in 
X saturates and one finds the fully nonperturbative result 

<?i(a;)«-l/4£, g2ix) ^ i-vofySs (31) 

and G^-*) - VX. 

In the case of spontaneous symmetry breaking at tree- 
level, m^g < 0, IR fluctuations restore the symmetry, 
resulting in a positive effective square mass cx A; see 
Eq. (fTTj). One has, assuming A <^ |m-dsl ^ 

e = Aeff/lmisI, TTp = 2\mls\ (32) 

and thus e — Hp/d 3> e. In this case, we find 

5i(x)=£e(x) and g2{x) =^ e{--qof'' Cl{x). (33) 

We see that loop effects are not suppressed by powers of 
the coupling. In the deep IR regime 7rp|a;| > 1, one has 

5i(x) «-!/£-, g2(x) « (-?7o)'7e" (34) 

which exhibits, again, IR enhancement: G'*' ~ -^/l^dsl- 
Finally, we note the specific momentum dependence of 
the 52 loop contributions in ([23l) . (f24| which is singular 
whenever the sum of any two momenta approaches zero, 
Kij — )■ 0. This loop contribution thus gives a distinct 



signature from the tree- level one and, in fact, provides 
the dominant contribution for such momentum configu- 
rations. This is a direct consequence of the IR behavior 
(HH) of the nonlocal four-point vertex p^ . In the present 
case, the latter is given by the two-point function ([13]) of 
the operator (jp' which, at large momentum separation is 
essentially that of a free scalar field of mass « d{£+£), 
as noticed in [l9| . Loop contributions can thus be seen 
as describing the exchange of a light (composite) scalar 
degree of freedom, whereas the tree-level contribution is 
a contact term. Similar contact and exchange diagrams 
appear in tree-level calculations of nouGaussian trispec- 
trum for primordial density perturbations in infiatio nary 
models and lead to specific nouGaussian signatures [23|- 

In conclusion, we have obtained an analytic expres- 
sion of the nouGaussian four-point correlator of an 0{N) 
scalar field in the large- limit, which shows that dS IR 
effects can lead to nonperturbative enhancement of loop 
contributions. We believe our result adds to the under- 
standing of the nontrivial IR physics of light scalar fields 
in dS space. Although the present model is too simplis- 
tic for realistic cosmology applications, our calculation 
demonstrates how IR effects may spoil the usual pertur- 
bative expectations. Possible implications for inflation- 
ary cosmology need to be investigated. 

We acknowledge interesting discussions with F. Gau- 
tier and R. Parentani. 
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